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ABSTRACT 

A simple algorithm is employed to deproject the two dimensional images of a pilot sample 
of 12 high-quality images of edge-on disk galaxies and to study their intrinsic 3 dimensional 
stellar distribution. We examine the radial profiles of the stars as a function of height above 
the plane and report a general trend within our sample of an increasing radial scalelength with 
height outside of the dustlane. This could be explained by the widespread presence of a thick 
disk component in these galaxies. In addition, the 3-dimensional view allows the study of the 
vertical distribution of the outer disk, beyond the break region, where we detect a significant 
increase in scalelength with vertical distance from the major axis for the truncated disks. 
This could be regarded as a weakening of the "truncation" with increasing distance from the 
plane. Furthermore, we conclude that the recently revised classification of the radial surface 
brightness profiles found for face-on galaxies is indeed independent of geometry. In particular, 
we find at least one example of each of the three main profile classes as defined in complete 
samples of intermediate to face-on galaxies: not- truncated, truncated and antitruncated. The 
position and surface brightness that mark the break location in the radial light distribution are 
found to be consistent with those of face-on galaxies. 

Key words: galaxies: photometry - galaxies: structure - galaxies: fundamental parameters 
- galaxies: evolution - galaxies: formation galaxies: individual: ESO 380-019, ESO404-018, 
FGC2339, IC4871, NGC522, NGC 1596, NGC3390, NGC4179, NGC5290, NGC5290, 
NGC5981.UGC 10459 



1 INTRODUCTION 

To study galaxy formation and evolution it is often helpful to imag- 
ine that galaxies are made up of a small number of separate com- 
ponents. Surface photometry of galaxies provides a very simple 
way of parameterising these various components. By modelling the 
surface-brightness distribution of a galaxy and thereby parametris- 
ing the individual components, we obtain a common ground to 
measure, compare, and sort large samples of galaxies. 

Owing to their flat geometry the random orientation of a spe- 
cific galaxy disk provides us with different projected views of 
the, presumably, common underlying 3 dimensional distribution. In 
particular, one gets complementary views from observing face-on 
or edge-on galaxies. In the face-on view, one studies the radial dis- 
tribution of the starlight discussing, e.g., non-axisymmetric features 
such as bars, rings, or spiral arms. In the edge-on view on the other 
hand, one is able to study the vertical distribution - e.g., disentangle 
thick disk, thin disk and stellar halos. Due to a longer line-of-sight 
projection, the edge-on view allows in principle the exploration of 
the radial profile to larger distances. 

Early on, from face-on galaxy observations, several authors 
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(e.g. |Patterson|[l94il Ide Vaucouleurslll959h noticed that the radial 
light distribution of stellar disks is well described by a simple expo- 
nential, a lbeit the lack of a clear physical explanation for this fact. 
Later on. iFreemanl i 1 970h called these disks Type I - as opposed to 
some other galaxies called Type II - showing a break in their radial 
profile with a downbending steeper outer region, or equivalently, a 
flatter inner region. 

Using the edge-on view, Ivan der Kruitl J 1979b discovered 
that the exponential decline does not continue to the last mea- 
sured point, but is truncated after several radial scalelengths. Af- 
terwards, with the advent of CCD technology, several groups 



jBarteldrees & Dettmarl 1 1994 IPohlen et alj boOObl; iFlorido et al.l 
1200 ll ; Ide Griis et al.ll200ll ; |Pohledl200ll : iKregel et alj|2002h stud- 
ied larger samples of edge-on galaxies to find out if and at what 
radius this truncation occurs. 

More recently. lErwin et alj ( 12005b . lErwin et alj fcOOfjb . and 
IPohlen & Truiillol i l2006h have analysed surface brightness profiles 
of face-on to intermediate inclined galaxies and revised our view 
on disk truncations. They identified three basic classes of surface 
brightness profiles depending on the apparent break features or lack 
of them: 1) The well known Type I that has an exponential profile, 
with no break. 2) Type II with a 'downbending break'. 3) A com- 
pletely new class, called Type III, also described by a broken ex- 
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ponential but with an upbending profile beyond the break. In this 
view th e classical truncatio ns d escribed initially in the pionee ring 
worko f l van der Kruiild 19791) and| van der Kruit & Searle ( 1981) are 
in fact part of a generalised Freeman Type II class (Freeman l 19701) . 
In this picture truncated galaxies show a surface brightness profile 
with a shallow inner and a steeper outer exponential region sepa- 
rated at a relatively well defined break radius. 

These r esults are in agreement with those reported by 
IPohlen et alj < l2002ah . who obtained deep photometry of three face- 
on galaxies and clearly demonstrated that the sharp cut-of fj] fit- 
ted to profiles of edge-on galaxies (e.g. by van d er Kruit & Searle 
1 198 it iBarteldrees & Dettmaril 1994k IPohlen et alf2 000b ) are in fact 
not co mplete but be t ter de scrib ed by a broken exponential. Even be- 
fore, |deGri]setal] feOOlh and lPohlenl ( 1200 lb noted that for edge- 
on galaxies the surface brightness distribution does not disappear 
asym ptotically ('vertically' ) into the background noise. Interest- 
ingly, iHunter & Elmegreenl d2006b find a similar broken exponen- 
tial shape for many objects in their large sample of irregular galax- 
ies. 

The redefinition of the truncation as a point at which the two 
exponential fits meet solves, among ot her things, the apparent con- 
tradict ion between simulations by iHotlev- Bockelmann & Mihos 
(2001) an d the edge-on obser vations of IBarteldrees & Dettmad 
1 1994b and lPohlen etai] {2000a). The latter authors find that 'trun- 
cations' are as 'close' as twice the radial scalelength from the centre 
in contradiction with the simulations. Furthermore, 'early' trunca- 
tions (i.e., at R ^ 2h) of edge-on galaxies have been always ques- 
tioned by the 'face-on viewers', who had never observed a face-on 
galaxy with a sharp truncation at a distance of twice the scalelength. 

Despite the line-of-sight enhancement of the surface bright- 
ness level in edge-on galaxies, star counts in nearby galaxies are 
able to provide significantly deeper profiles with much higher S/N . 
Using this completely independent method. [Ferguson et al.1 ( 2006) 
followed the profile of the face-on galaxy M 33 down to /iii m ~ 30 
Tmag/n" , which is currently unreachable with surface photom- 
etry. Their result shows that the broken exponential description is 
valid down to the last bin above the noise level and an actual edge 
of the disk is not found, at least for this galaxy. 

Nevertheless, there still seem to be some open issues for a 
complete reconciliation o f the renewed face-on with the 'old' edge- 
on view dde Zeeuwll2007h as to provide a consistent picture of the 
three-dimensional structure, for example: 

• Is the broken exponential structure, now well observed for 
face-on galaxies, quantitatively consistent with the truncated pro- 
files observed for edge-on galaxies? 

• Do we see a different kind of truncation in the deeper views 
provided by the edge-on galaxies? 

• Can we identify the edge-on counterparts of the three 
main types of surface brightness profiles (Type I, Type II, and 
Type III) observed for fa ce-on galaxies bv lErwin et al.1 d2006b and 
IPohlen & Truiillol d2006b ? 

• Does the frequency of different profile classes found for the 
face-on galaxies match the edge-on results? 

• Why did the former edge-on surveys, with large and deep sam- 
ples, not discover the new class of Type III profiles? 

To answer these questions and disentangle the crucial problems 



1 It is interesting to note that in this concept the truncation radius is often 
equal to the last measured point and therefore much further out compared 
to the break radius of a broken exponential fit. 



with the line-of-sight integration (see e.g. IPohlen et al]|2004ab wc 
decided to deproject a sample of edge-on galaxies. This can only 
be done using some geometric assumption. Here we choose the as- 
sumption of axial symmetry. This will allow us to gain insight into 
the three-dimensional galactic stellar disk structure by bringing to- 
gether the classical face-on and edge-on views. 

The vertical structure may well be the key to understand the 
nature of the different breaks in the surface brightness profiles, 
by setting additional constraints on the various physical explana- 
tio ns put forward. For example, some recent star formation mod- 
els jElmegreen & Hunter 2006; Li et al.l2006b are now able to pro- 
duce the observed bro ken exponential structure, but, alternatively, 
iDebattista et al.1 |2006) also find downbending breaks using purely 
collisionless iV-body simulations. Once the nature of the galactic 
stellar disk is determined it will be fully justified to utilise these 
breaks for c o mparison of galaxies at v ari ous redshifts as done b y 
IPerezI d2004b , lTruiillo & PohlerJ d2005b or lTamm & Teniesl d2006b . 
an d explain their appearan ce in cosmological simulations as shown 
bv lGovernato et al.1 12007b . 

The remainder of this paper is organised as follows. Section|2] 
describes our sample selection and data. In Section|3]we introduce 
our method to deproject the edge-on galaxies. The results are given 
in Section [4] we discuss them in Section [5] and conclude in Sec- 
tion^ In Appendix lAlwe give detailed comments for all our galax- 
ies and show in Appendix [B] radial surface brightness profiles and 
isophote maps for all galaxies. 



2 DATA 

All optical images used in this paper were taken as part o f a PhD 
study on the radial structure of galactic stellar disks by IPohlenl 
d200lh . Originally, the disk galaxies (-2 < T < 7) were cho- 
sen according to the allocated observing time, the observatory 
(north/south), and the available FOV (D 2 s ^2") to meet the fol- 
lowing morphological selection criteria verified by using images 
from the Digitized Sky Survey (DSS): 1) The galaxies are edge- 
on with i <; 86°. 2) Undisturbed with smooth photometric features 
(see below). 3) They look similar to some prototypical cases like 
NGC4565 or IC2531, this is done to make it possible to consis- 
tently fit a simple disk model. 

Galaxies with the following characteristics were rejected: 1) 
Apparent spiral arms, indicating a lower inclination. 2) A signifi- 
cantly asymmetric or disturbed disk, indicating strong interaction. 
3) Two-sided or significantly one-sided warped disks. 4) Galax- 
ies apparently dominated by the light of their bulge component 
(i.e., there are almost no Sa/Sab galaxies, but several lenticulars). 
5) Galaxies that showed only a faint, patchy, not well defin ed disk. 
For m ore details of the original sample we refer the reader to lPohlenl 
d200lh or lPohlen et all d2004bl) . 

The present sample of 11 galaxies is drawn from this base 
sample of 72 galaxies as a pilot sample covering a wide range of 
Hubble types and apparent vertical distributions (thin and thick) 
while trying to pick galaxies with very few disturbing fore- or back- 
ground objects. We add an additional NIR image of one of the 
galaxies to observationally assess the influence of the dust. Global 
properties of the finally chosen galaxies are given in TableQ] 

The images were obtained in six different observing runs be- 
tween 1998 and 2004. Two runs (El and E2) at the Danish 1.54m 
telescope and one (E3) at the 3.6m NTT telescope, both at the Eu- 
ropean Southern Observatory in Chile. One run (CI) at the 1.23m 
telescope and one (C2) at the 3.5m telescope both at Calar Alto in 
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Table 1. The 1 1 galaxies of our pilot survey. 

(1) Principal name, (2) right ascension, (3) declination, (4) RC3 coded Hubble-type, and the (5) Hubble parameter T (all taken from de Vaucouleurs et al. 1991), 
(6) diameter in arcminutes, (7) heliocentric radial velocities, (10) maximum deprojected rotation velocity of the gas, and the B-Band absolute magnitude (11) 
(taken from LEDA). According to the heliocentric radial velocities corrected for the Local Group infall into the Virgo cluste r (8) from LEDA, we e stimated the 
(9) distances following the Hubble relation with the Hubble constant from the HST key project of Hq = 72 km s~ x Mpc _1 (Freedma n et alj200ll) . 



Spain. And finally one run (LI) at the Hall 42" (1.06m) telescope 
of the Lowell Observatory at Anderson Mesa (USA). 

During the two runs with the 1.54m Danish telescope it was 
equipped with the DFOSC camera and the C1W7/CCD which is a 
2kx2k LORAL chip providing a field size of ~ 13' and a scale 
of 0.39"pixer\ At the NTT we used EMMI, equipped with a 
2kx2k Tektronix chip providing a field size of ~ 9' and a scale 
of 0.27" pixel -1 . The run at the Calar Alto 1.23 m telescope was 
done in service mode with the Site#18b chip, a 2k x 2k SITE CCD 
with 24 /im pixel size, providing an unvignetted field of ~ 10' and a 
scale of 0.5"pixel -1 . At the 3.5m telescope we used the wide-field 
near infrared camera OmegaPrime with a lkx lk pixel HAWAII 
HgCdTe array by Rockwell, providing a field of ~ 7' at a scale of 
0.4"pixel _1 . During the run at Anderson Mesa the 42" telescope 
was equipped with the NSF CCD used with 2x2 binning providing 
a field size of ~4.9' and a scale of ~ O^S'^ixeP 1 . For the optical 
imaging the available standard Bessel R and V filter at ESO, and 
Johnson R at Calar Alto and Anderson Mesa are used. The NIR 
image with OmegaPrime is in the K'-band. Table[2]summarises the 
detailed observational parameters for each image. 

For more details about the data reduction (i.e. flatfielding, mo- 
sa icing, and sky subtraction) and phot ometric calibration we refer 
to lPohlenl j200ll) : |Pohlen et al.l J2004bl). De t ails of the NIR data re- 
duction can be found in lvan den Berg et al] ( l2006l) . 

2.1 Preparation of the images 

Before deprojecting our galaxy images we first have to prepare the 
images, here we describe the procedure followed in this work (see 

Fig.[Q- 

As a first step we mask all the foreground stars and the back- 
ground galaxies. To do this semi-automatically, we use SExtractor 
(Berti n & Arnoutsll 19961) to create an object catalogue. Depending 
on the measured ellipticity and the provided star/galaxy classifier 
we separated stars from galaxies, which in turn get assigned circu- 
lar masks with radii depending on their measured flux. The output 
is a DS9 region file which is overlayed on the galaxy image and 
visually inspected for each of the masked objects. For highly el- 
liptical background galaxies the mask is changed to an ellipse. In 
addition we have to adapt the radii of the masks. Typically one opts 



for a conservative masking (i.e., prefer to mask more than necessary 
and leave 'larger holes') to make sure that all non-galaxy features 
are removed. 

However, since our deprojection method cannot deal with 
masked areas (actual holes) we have to interpolate the flux across 
our mask. By doing this, one is forced to chose the mask radius 
as small as possible in order to not introduce too much unwanted 
structure. This approach clearly leaves residuals on the masked, in- 
terpolated image which will be discussed later (cf. Section [3~7l >. 

In a second step the galaxy image is rotated to align the galac- 
tic plane parallel to the major axis using the smallest angle of rota- 
tion according to their true position on the sky (chip). The rotation 
and centering is done by blinking the adjacent sides and minimis- 
ing the difference by eye. Depending on the intrinsic symmetry this 
can be done down to an accuracy of ±0.5° . 

The third step is to cut the galaxy into 4-quadrants and av- 
erage them. This procedure increases the S/N and minimises the 
residuals from the interpolated, masked regions and gives a more 
regular averaged galaxy. The influence of the dustlane and intrinsic 
asymmetries will be discussed later (cf. Section [3~7l (. 

The final step is smoothing the image. Since our deprojection 
code only handles one radial cut at a time, we choose the kernel 
size for the Gaussian smoothing to match the vertical increment 
between successive deprojections in order to get fully independent 
cuts but using all available rows. 



3 METHOD 

In this section we describe and test the method used to recover the 
3D structure of the available disk galaxies. We start with describing 
the deprojection algorithm (Section [3~Tt and discuss the noise prop- 
erties of the reconstructed 3D distribution, especially after the orig- 
inal image has been smoothed with a Gaussian kernel (Section [3~2l . 
To further reduce the noise we applied an additional denoising 
(Section [3.3| >. In Section [3^41 we discuss in detail our choice of fit- 
ting function for the radial surface brightness distribution. Then we 
describe how the radial profiles of the recovered distribution for a 
given z are actually fitted (Section [3.5t . In Section [3~6l we apply 
the method and the fitting procedure on a large sample of artificial 
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Table 2. Observing log for the individually combined images. 

(I) Galaxy name, (2) filter, (3) observing date, (4) site, (5) the total coadded on-source exposure time, (6) seeing, as measured on the final combined image, (7) the 
number of individual images with their individual exposure times. 
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Figure 1. Preparation of the images. 

Contour maps of NGC 598 1 showing the four main steps in image preparation. The upper left panel shows the rotated, original image. The upper right panel 
shows the image after masking of fore-and background sources and interpolating across the masks. The lower left panel shows the galaxy after averaging the 
four quadrants. The lower right panel shows the galaxy after Gaussian smoothing. The equidistant contour levels plotted are in all cases the same, except for 
the smoothed version where two additional outer contours are added. 



galaxy images in order to test the method's ability to recover the 
underlying 3D distribution. In this set of artificial images the influ- 
ence of the inclination of the galaxy and the existence of dust on 
the recovery of the 3D distribution is tested. Finally, we estimate 
the effects of our image preparation on the resulting profiles using 
real galaxies in Section [3~71 



3.1 Deprojection 

The deprojection of axial ly symmetric quantities is a classical prob- 
lem in astronomy (e.g., iLucvl 119741) . It is considered for deter- 
mining the 3D ste ll ar orb it struc tu re in elliptical galaxies (e.g., 
Dehnen & Gerhard 1 19931 1 1994 iBinnev, Davies & Illingworth 
19901) . and in determining clus ter 3D structure from X-ray ob- 
servations dFabricant et~alj|l984l) and more generally from X-ray, 
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Sunyaev-Ze ldovich and cluster weak-lensing data (Zarou bTet al.l 
1998Ll200ll) . In the case of edge-on galaxies. iFlorido et al.1 j200lL 



2006) have applied an inversion method to obtain deprojected pro- 



files before, but without studying a dependence with the vertical 
distance from the plane. Here we apply a simple analytical method 
that employs the, so called, Fourier Slice Theorem in order to re- 
cover the intrinsic 3D stellar distribution from our sample of high 
quality 2D edge-on disk galaxy images. 

The disk galaxy images are deprojected as follows. We adopt 
the convention that bold-face symbols denote 3D quantities (e.g., 
k = (k x ,k y ,kz)). Let the observer's coordinate system be de- 
fined with the Cartesian axes (x, y, z), with the y axis aligned with 
the line-of-sight. The disk galaxy axis of symmetry lies along the 
z-axis which is perpendicular to the line-of-sight. Let I(x, z) = 
J A(x) dy denote a projected quantity (image) of the source func- 
tion A. The 3D Fourier transform of the source function is related 
to the image by 



\(k x ,0,k z ) = 



J-i(.k x x+k z z)] 



I(x, z) dx dz 



(1) 



Note that functions marked with ~ are in Fourier space. 

This relation could be easily understood since the galaxy has 
a rotational symmetry, i.e., A(r) = A(r, z) and its Fourier trans- 
form is A(k) = A(fc, k z ) where k — ^Jk x + ky. The deprojection 
is then obtained by inverse Fourier transform of I(k x , k z ) and re- 
placing k x with k, namely by using the following equation: 



X(r,z) 



(2tt)2 



k ^I{k,k z )J (kr)kdkdk z 



(2) 



Here Jo is the zero order Bessel function. 

Notice that for a specific radial slice, i.e., a cut at a given zo 
and parallel to the major axis, the integral I0 can be written simply 
as, 



A(r, 20) 



1 

W) 



T(k, zo)Jo(k r)kdk . 



(3) 



WhereT(fc, zq) is the Fourier transform over the image x-direction. 
This way of writing the equation enables us to easily deproject cer- 
tain radial slices of the image, or collapsing the whole z-direction 
of the image to recreate the equivalent of a face-on view of the 
galaxy. 

In principle, it is easy to obtain an equation similar to (O for 
the general case of a disk with an arbitrary inclination angle. How- 
ever, one can show that for such cases there is a cone in k-space 
inside which there is no information on the 3D structure in the 
image. The general case is beyond the scope of this paper and is 
not needed for the range of inclinations we restrict our selves to 
i.e., to i <; 86° (see section|3.6> . The interested reader is referred to 
IZaroubietalUl 99811200 lh . 

Furthermore, one can show that the Fourier Slice Theorem 
approach for inverting edge-on galax ies is actually another wa y 
of writing the inverse Abel transform iBinnev & Tremaine"lll987l) . 
The advantage of the Fourier space language here is that it avoids 
taking direct derivatives of noisy data, hence, making the inversion 
more stable. 

The applied Fast Fourier Transform (FFT) assumes periodic 
boundary conditions which has to be satisfied by our images. 
Therefore, one has to pad the image with at least twice the number 
of the image pixels along each axis. However, since FFT of sharp 
transitions will give a strong aliasing signal in the deprojected slice, 
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Figure 2. Test of deprojection method for NGC 598 1 . 
In the upper frame the radial profiles of the symmetrised image (dashed- 
lines) and those of the deprojected-reprojected image (solid-lines) are 
shown. The uppermost profile is the vertically integrated profile, while each 
of the others represents a cut at increasing vertical distance from the ma- 
jor axis. The lower frame shows the percentage error one obtains from this 
procedure. 



we do not apply zero padding. Rather we padded the image with 
white noise background with the same rms as the noise in the ob- 
served image. The residual small differences at the edge changing 
from the artificial white noise to the actual noise in the background 
of the image is washed out by the final wavelet denoising and Gaus- 
sian smoothing of the image (see Section [3~2l and [3~3l . Further tests 
with more careful tapering of the image's edge with a Fermi-Dirac 
function without adding noise showed that there is no significant 
influence on the deprojected profiles. 

As a zero-order test, we deproject the image of NGC 5981 and 
then reproject it back. The result is that we obtain profiles that are 
identical to those of the original image (cf. Fig.[2} down to < 0.5% 
error level. Most of this error actually comes from the reprojection 
- not the deprojection - procedure. The reason is simply that the de- 
projection is a differential quantity that depends on the intensity of 
the image at a given radius, whereas the reprojection is an integra- 
tion procedure that will depend on the intensity of galaxy at every 
radius including radii beyond those explored in the image (i.e., af- 
ter the noise dominates over the signal.) A more thorough testing 
of the method in which we compare the deprojected 3D structure 
with the original one will be presented in Section [3~6l 



3.2 Noise properties and smoothing 

To explore the noise properties, let us suppose that the image only 
contains white noise, with rms a, i.e., (I(xi, zi)I(x2, Z2)) = 
(2tv) 2 S(x 1 - x 2 )8(zi — Z2) a 2 . The noise correlation function 
in the deprojected image is given by 
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(A(ri,zi)A(r 2 ,2a)) 



22 



(27T) 



fc 2 J (fcri)J (fcr 2 )dfc.(4) 



Not surprisingly, the high frequency components dominate and 
hence some smoothing must be employed. The analytic nature of 
the correlation matrix is an attractive feature of the method. 

Another nice feature of the method is that any isotropic 
smoothing function applied to the image will depend only on the 
modulus of the radial wavevector, jkj, in the deprojection. This can 
easily be seen from the previous equations. A simple and practi- 
cal Gaussian smoothing filter, with smoothing scale R s , gives an 
auto-correlation of the noise of the form 



M{r) = (\{r,zf) 



(5) 



where 2F2 is a hyper- geometric function. This has a particularly 
simple asymptotic form, 



lim Af(r) 



(6) 



Since we are interested in the large scale features of the recon- 
struction a Gaussian smoothing is applied on the image. The choice 
of the FWHM of the smoothing kernel is discussed in l2.ll 

Despite these simple analytic features of the noise properties, 
we sample the noise properties of the smoothed 3D structure cre- 
ated by the method with deprojecting noisy slices of the image 
(with a given z well outside the galaxy). For each image we de- 
project some 10 randomly selected lines, average their rms and use 
that as our noise estimate. The profile of the reconstructed noise 
rms is shown in every set of reconstructed profiles (see Fig.[5}- 

The noise is dominated not by pixel to pixel variations, but 
by any residual large scale asymmetries in the background as well 
as residuals from the masking (cf. Section |2~TV The latter cause 
the somewhat bumpy structure of the resulting deprojected profiles 
(see Fig.[5]l. 

The seeing acts in effect as an additional Gaussian smoothing 
of the image. This type of smoothing flattens the brightness profile 
around the center but has very little influence on the scalelength 
measurements. 



galaxies ( Hunt er & ElmegreenI l2006h suggest that the broken ex- 
ponential is a good fitting function for the radial light distribution 
(cf. SectionQJ. Technically, this results in a discontinuity between 
the inne r and outer slope, which is unphysical and of course not 
implied. IPohlen & TruUliol d2006h find an extended transition or 
break region, whic h is not as sharp as measured before (e.g., by 
IPohlen et alj2002d) . From their sample only a few clean and undis- 
turbed truncated cases (e.g., NGC5300 or NGC7437) exhibit a 
very sharp transition between the inner and outer disk region. Most 
often their Type II breaks appear with a more gradual transition 
zone of a size up to ~ 4 kpc. Despite the fact that the exponential 
function in itself has no clear physical explanation, the broken ex- 
ponential is a reasonable fitting function that describes the shape of 
the profiles in a simple and uniform way. 

The key point of this description is that the outer disks of 
truncated galaxies - regardless of any transition zone and the ex- 
act shape of the outer profile - makes up a significant fraction of 
the physical size of the their stellar disks. For example, the ratio 
between the inner and t he outer disk is typically 2 to 1 in radius 
jPohlen & Truiilloll2006h . therefore in terms of area the outer disk 
is even slightly l arger. In some cases authors ignore the outer region 
altogether (e.g. Jvan der Kruit & Searldl 198 it IPohlen et al.l l2000b) 
and use an exponential with a sharp cut-off as their fitting function 
of choice. However, this is not the case for edge-on galaxies in gen- 
eral as clearly appare nt from the examination of the radial profiles 
of, e.g., NGC 5907 in Ivan der Kruit & Searld ( fl98lh or IC 4393 in 
|Pohlenetal] j2000b). We see the same for several galaxies used in 
our study. For example, the major axis profiles of ESO 380-019, 
ESO404-018, NGC 5981, UGC 10459, and NGC 522 (see Fig.lBTl 
in Appendix |Bj a re obviously better describ ed by a broken expo- 
nential fit (Pohlen 200 H; IPohlen et alj|2002bh. Many further exam- 
ples in the literature (e.g. , see IJensen & Thuanll 1982l;ISasakill 19871: 
iNaslund & Jorsated[i997l ; lBvunlll998l ; |Pohlenlf200l| - |de Griis etail 
2001) show that the broken exponential shape is indeed a very good 
fitting function even for edge -on galaxies. However, w e should 
mention that some authors (e.g. lFlorido et al .l200lll2006h . still pre- 
fer the description of a smooth, but complete truncation. 

To summarise, since the profiles of face-on galaxies are well 
parameterised by a broken exponential model, and those of edge-on 
galaxies seem to be, it is useful to apply the same parametrisation 
for the profiles of the edge-on galaxies studied here. 



3.3 Wavelet Denoising 

In order to reduce the apparent bumpy structures along the radial 
profiles we apply an extra denoising algorithm based on wavelet 
expansion. The employment of this algorithm is not necessary and 
we tested that it does not change the quantitative results of the de- 
projection in any significant way, it however reduces the apparent 
"bumpiness" in the radial profiles. 

The denoising schem e adopted here follows the method pro- 
posed by iDonohol dl992l) (see also [Donoho & Johnstond 1 19941 ; 
IZaroubi & Goelmanl 120000 where the wavelet coefficients are re- 
duce with a constant equal to the rms of the white noise. This proce- 
dure is known as the soft thresholding of wavelet coefficients. The 
wavelet functions used here are the Daubechies' discrete wavelets 
dDaubechieslll988h with 12 parameters (DAUB 12). 



3.4 Profile shape and fitting function 

The recent studie s of profile shapes for face-o n disk galaxies 
jErwin et alj |2005| , 120061 : IPohlen & Truiillol I2006T) and irregular 



3.5 Parameter fitting 

In order to quantify the behaviour of the recovered 3D stellar dis- 
tribution of the galaxies, we fit the deprojected radial profiles, for a 
given z, with a simple broken exponential function of the form 



ft{R) 



Mo,m + 1-086 * R/fn 
Mo.out + 1.086 * R/h, 



if R ^Rbr j 
, if J? >Rbr ■ 



In principle, the model includes 5 free parameters, the in- 
ner and outer scalelengths h[ n and h out , the inner and outer cen- 
tral surface brightnesses, fio,m and /Lto.out and the location of the 
break radius, 7?br , which marks the boundary between the inner 
and outer regions. How ever by using the downhill simplex method 
dNelder & Meadl 965) with initial starting values of only the scale- 
length and central brightness parameters we avoid using the break 
radius as a free parameter. This is achieved by deciding a-priory 
that the inner slope is either shallower or steeper than the outer 
one. This enables the break radius to be determined by the crossing 
point of two exponential fits found by the algorithm. We find that in 
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all cases, except NGC 1596, the inner profile is shallower than the 
outer profile. In addition, when the ratio between the two fitted ra- 
dial scalelengths is about 1 ± 0.2, we fit it with a single exponential 
profile of the form fi(R) = no + 1.086 * R/h. 

In order to avoid including the bulge and/or bar regions in the 
fitting procedure, i.e., avoid a detailed bulge to disk decomposi- 
tion, we used the "marking the dis k" method with fixed inn er, bi, 
and outer, &2, boundaries (see e.g. IPohlen & Troiilldl2006f) . The 
inner boundary is chosen to exclude the region that is obviously 
dominated by the bulge plus bar (shoulder) component. The outer 
boundary is set at the point at which the surface brightness profile 
is comparable to the noise level (see Section [3T2l , In the case of 
an extended bump clearly visible in the profile, caused for exam- 
ple by a residual from the masking, the boundary is moved further 
inside. Sometimes a sharp rise or drop appears beyond the outer 
boundary of the recovered profiles (see, e.g., the face-on equiva- 
lent profiles of IC4871 or NGC 5290 R-band in Fig. [5). This is 
caused by the Fast Fourier transform (FFT) operation used in the 
method. The FFT procedure requires periodic boundary conditions. 
To avoid contamination of the profile at hand by these conditions 
we increase the size of the image by a certain factor. However, this 
periodicity still shows up in the noise dominated regions and some- 
times close enough to, but never within, the region of interest. 

The five starting values for the downhill simplex minimisation 
are obtained by a single exponential fit to the inner or outer region 
with an approximate break radius for each galaxy. 

3.6 Tests with artificial data 

3.6.1 Model grid 

To study the influence of the dust and possible non-90° inclina- 
tion on the scalelength measurements in our deprojected profiles 
we used a grid of 96 mock images (400x400 pixels) kindly pro- 
vided by Simone Bianchi. Using a Monte Car lo radiative transfer 
method iBianch i! Ferrara, & Giovanardilll996h . he calculated sin- 
gle disk models at a range of edge-on inclinations (80° — 90°) with 
different input parameters for the stellar and dust dis tribution. A 
similar approach is adopted bv lMollenhoff et alj |2006) to measure 
the influence of dust and inclination for less inclined galaxies. 

The chosen inclination angles are 80°, 85°, 86°, 87.5°, 89°, 
and 90°. The stellar disk is a 3D single disk model with expo- 
nential profiles both along the radial and vertical direction. The 
radial scalelength is h — 33.4 pix, and the vertical scaleheight is 
zo = 0.125ft. The only parameter of the stellar distribution that is 
varied is the position of an outer (sharp) truncation of the light dis- 
tribution. The truncation parameter has 4 values, 10ft, 5.4/i, 3Ah 
and lAh. The first value mimics an infinite exponential profile, and 
the other three values represent a long truncated, an intermediate 
truncated, and a short truncated stellar disk. 

For the dust distribution a set of four different cases with in- 
creasing influence on the measured R-band light distribution are 
adopted. The parameters for t hese cases are within the range mea- 
sured bv lXilouris etail dl999f) . The first case is a dust free (trans- 
parent) model which we use as our reference case. The second 
case, which we call the "best case", has the following dust pa- 
rameters: tr = 0.20, ha/h — 1.0, and Zd/zo = 0.30. Here 
tr is the R-band optical depth, hd is the dust scalelength, Zd is 
the dust vertical scaleheight and zo is the stellar scaleheight. The 
parameters of the third case, called the "intermediate case", are 
tr = 0.54, ha/h* — 1.4, and Zd/z, =0.55. Finally, the parameters 
of the fourth case, which we call the "worst case" are tr = 0.65, 



hd/h t = 1.5, and Zd/ z* — 0.8. For the sake of simplicity the dust 
disk is truncated at the same radius as that of the stars. In addition, 
we have tried several other cases with variable dust truncations, 
however those did not show any significant difference with respect 
to the ones presented here. 

The number of images used to test the dust and inclination 
influence on the reconstruction is then doubled (except in the 90° 
cases) by using two different methods of symmetrising the far and 
near side of the galaxy, insofar as the dustlanes are concerned. In 
the first case we average all four quadrants, whereas in the second 
we only average over the two far side quadrants which have less 
dust influence and ignore the two near side quadrants with the dust 
lane. 

The artificial images are 3x3 pixels carbox smoothed. The de- 
projection is carried out on the vertically integrated profile, here- 
after, the face-on-equivalent view, as well as on the major axis and 
two higher cuts at z = 1.7zo and z — 3.4zo- Since the input profile 
in this case is sharply truncated, we apply only simple exponential 
radial fits ranging from the center to 3 times the radial scalelength, 
except for the short truncated disk where we go out to one scale- 
length. From these fits we obtain the scalelength and central surface 
brightness of the deprojected 3D distribution. 

3.6.2 Influence on the measured parameters 

To visualise the influence of the various dust distributions we plot- 
ted the relative deviation in the recovered scalelength as a function 
of the inclination angle (Fig. [3j. Here the relative deviation is de- 
fined explicitly as 100 X (/itrue— /lrccovcrcd)/frtruc- 

From Fig. [3] it is obvious that using only the two dust-free 
quadrants, shown as dashed-crossed lines, improves the quality of 
the fit. However, above 88° the gain is only minimal and coun- 
tered by the advantages of an increased S/N - especially important 
to reach large vertical distances - and by minimising the residu- 
als from masked objects. Therefore, for the image preparation of 
the real data we average over all four quadrants, even in the clear 
presence of dustlanes. 

As expected our deprojected, face-on-equivalent, scalelengths 
are very close to the true, 3D input scalelength in the absence of 
dust. There is almost no influence of the non edge-on geometry 
(down to 80° ) of the actual model galaxy which we treat in our de- 
projection to be always at 90° . For the model galaxies with dust, 
we see that the error on the scalelength is almost always systemati- 
cally smaller for galaxies with increasing deviation from the perfect 
edge-on geometry. 

We also notice, from the dust-free case, that for the short- 
truncated disk the deprojection improves the fits quite noticeably, 
since it disentangles the truncated part of the profile and avoids fit- 
ting for it, as is the case with the projected edge-on view. For the 
latter we know that they result, slightly depending on the region 
used for fitting, in a syste matic overestimation of the scalelength 
in the range of 20% (e.g., |Pohlerj|200lLlKregel et all2002klKregeil 
l2003h . 

However, we note that our face-on-equivalent scalelengths are 
still systematically overestimated by typically, 15% — 30% due to 
neglecting the existence of the dust component. Naturally the ef- 
fect is far the worse on the individually deprojected major axes 
(Ziin overestimated by more than 50%), which probably make only 
sense for dust-deficient SO galaxies or in the dust-reduced NIR 
bands. However, the narrowness of the dustlane works in our ad- 
vantage, so that above 1.4 times the vertical scaleheight the error 
in the measured scalelength already drops to less than 10% (for 
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Figure 3. The systematic deviation of the measured scalelength with inclination and dust influence. 

Plotted is the deviation of the measured from the known input scalelength (in percent) from fitting the 96 mock images with the same setup compared to the 
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Figure 4. Influence of asymmetries: Averaging quadrants 

The deprojected surface brightness profiles of NGC 5981 using four different ways to combine the four quadrants. The upper left and upper right panel show 
the result using only the two quadrants of the south-east and north-west side respectively. The lower left panel shows the profile using only a mean of the 
two dust-free quadrants. The lower right panel finally shows the resulting profiles after averaging all four quadrants. In each panel, the dashed line (the upper 
profile) is the face-on-equivalent, radial profile of each galaxy. They are obtained by summing the light of the edge-on galaxy along the columns and deproject 
the result. The following dotted-dashed lines are the deprojected profiles of individual cuts from top to bottom starting from the major axis towards cuts at 
larger distance from the plane. The proper distance z from the plane is given in the caption (together with name and filter) of each panel in the upper right 
corner. The dotted line (the lower profile) is the deprojected profile of a noise (or source free) cut outside of the galaxy. Overplotted are on each profile two 
squares marking the fitting region for the disk model. The result of the broken exponential fits are overplotted on each profile as solid lines where the cross 
marks the break radius. 



i <; 87°). This means the scalelength measured on the individually 
deprojected higher vertical cuts are much better. 

The discussion of the results that we obtain from the observed 
galaxies takes the effects of the dust and inclination, that we learned 
from this section, into account. 



3.7 Influence of intrinsic asymmetries 

Dustlanes: In the last section we have estimated the influence of the 
dust on the measured flattening of the inner scalelength using our 
artificial images. Aside from a solely quantitative correction to be 
made to the measured scalelength in this case, one could argue that 
the dust in addition might be responsible for the observed broken 
exponential structure. In this picture the break would correspond to 
the point at which the dust lane becomes less important. However, 
in almost all cases the location of the break does not coincide with 
the, admittedly difficult to estimate, location of the visible end of 
the dustlane on the image. An obvious exception is the major axis 
profile (z = 0") of NGC 3390 where the dustlane extends out to 
R ~ ±55" and indeed causes an apparent break feature at this point 
(see Fig.O, hence it is excluded from the fit. In addition, we still 
find profiles showing a broken exponential structure at more than 
two times the stellar scaleheight above the disk, where the effect 



of the dustlane is minimal. Also NGC 4 179 does show a break, but 
has no dustlane at all. 

Furthermore, if one assumes that the outer slope is the intrin- 
sic, dustfree slope, the extrapolation into the center will in almost 
all cases come up too high (much brighter than the bulge peak). 
This would require far too much dust in these galaxies. There- 
fore, the dust quantitatively affects the measured inner scalelengths, 
however, is not responsible for the break in the profiles (cf. Fig.[4]l. 

Intrinsic small scale asymmetries: For real galaxies the residuals 
from masked objects will certainly affect the recovered profile. The 
problem is ostensibly aggravated by the averaging process of the 
four quadrants (cf. Section l2Tt , which will mirror any residual on 
one, to all other sides. However, one has to keep in mind that our de- 
projection method cannot deal with 'holes' , so we have to somehow 
fill these masked regions. For example, in the case of NGC 5981 
the extended companion superimposed on the upper left side of 
the bulge could be well traced in the contour lines of the averaged 
galaxy image (see Fig. []}. Under these worst case circumstances 
we made sure that this area is excluded from our fitting region. In 
the latter case we did not fit beyond z > 20". 

In addition, the final position of the break radius is strongly 
affected by the small scale asymmetries (star-forming regions, dust 
concentrations, residuals from masked objects) which vary in each 
of the vertical profiles. Whereas this only slightly affects the mea- 
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sured scalelength, the meeting point of the two exponential fits 
(break point) is more susceptible to such variations. 
Large scale asymmetries: We do not expect significant large scale 
asymmetries in our galaxies, since the input sample is morpholog- 
ically selected to be symmetric. NGC5981, which is asymmetric 
relative to the rest of the sample, is a good test case to estimate 
the effect of averaging quadrants from both sides. The asymme- 
try is clearly seen in the radial profiles (cf. Fig. IB 1 I t. The outer 
slope of the broken exponential fit is significantly shallower on 
the NW(right)-side, where also the contour map is more extended. 
Fitting only a mean quadrant obtained from the NW(right)- or 
SE(left)-side of the galaxy separately (cf. Fig. |4j we find differ- 
ences of 5% — 10% for the scalelengths of the inner region. In 
the outer disk (beyond the break) the differences are clearly larger, 
e.g., 40% on the major axis, due to the stronger asymmetry in the 
outer disk. However, both sides are equally well fit with the broken 
exponential function, so we do not introduce a bias on the validity 
of the applied fitting function by averaging the two sides. 



4 RESULTS 

Now we move on to present the results obtained by applying the 
deprojection method on our observed edge-on galaxies. Again here 
we first apply the method on the vertically integrated image. This 
results in a profile that is equivalent to the one obtained from az- 
imuthal averaging, for example by fits with concentric ellipses, of 
face-on galaxies. A direct comparison between the two views (see 
Section [3~0 is carried out in Section l4~T1 

Furthermore, we deproject a number of vertical slices to ex- 
plore the existence and properties of, for example, a thick disk com- 
ponent. Here the following questions are addressed: does the inner 
slope flattens with increasing z, as expected for some thick disk 
models? Does the break radius, if present, vary significantly with 
z? Finally, how does the vertical distribution of the outer disk, be- 
yond the break point, changes with z? These question are addressed 
in Section |4~2l 

The radial profiles obtained from deprojecting all the galaxies 
in the sample, together with the best fit broken exponential models 
are shown in Fig. [5] The results for the face-on-equivalent profile 
of each galaxy in the sample are tabulated in Table[3] 

4.1 The face-on-equivalent view 

4.1.1 Face-on vs. Edge-on: profile types. 

As mentioned earlier there are three basic classes of surface bright- 
ness profiles that hav e been found for face-on galaxies, Type I, 
Type II, and Type III dErwin et alj|2005l , 120061 ; IPohlen & Truiillol 
2006). Here, in the deprojected face-on-equivalent view, we also 
find at least one example for each of the three types of surface 
brightness profiles. 

All eight Sbc and later-type galaxies (ESO 380-019, ES0404- 
018, FGC2339, IC4871, NGC522, NGC5290, NGC5981 and 
UGC 10459 ) show clear Typ e II, downbending profiles. In compar- 
ison with lErwin et al. (2006) and lPohlen & Truiillol dixirl we can- 
not increase the level of classification (i.e., assigning sub classes) 
in order to explore the possible physical origin of the break, for ex- 
ample, star-formation threshold or bar related break. This is mainly 
because inferring accurate bar sizes in edge-on galaxies is almost 
impossible. 

It is interesting to note that the NIR K'-band profile of 



NGC 5290, typically assumed to be dust free, shows exactly the 
same behaviour as the R-band. Therefore, the dust and recent star 
formation, clearly more prominent in the R-band, do not play a role 
in the basic classification scheme applied here. 

The lenticular galaxy, NGC 4179, also exhibits a Type II pro- 
file with a ratio of inner to outer scalelength of hi n /h ou t = 
1.4. This ratio is somewhat lo w in comparison to the mean of 
{hin/h ou t) = 2.1 reported bv lPohlen & Truiillol d2006h for late- 
types, or the median of 2.4 reported by lErwin et al.l d2006l) for 
early-type galaxies. 

One galaxy, NGC 3390, although classified as an Sb galaxy, 
exhibits a dominant outer spheroidal component closely resem- 
bling an SO galaxy (see Fig. IB It . Its face-on-equivalent profile is 
consistent with being a classical unbroken profile, Type I. Never- 
theless, we can also fit to the recovered profile of this galaxy a 
broken exponential, a Type II. The scalelength ratio in this case is 
only hin/hout ~l-2, and so close to unity that in case of a face-on 
galaxy this would not qualify as a Type II. 

The profile of the second SO ga laxy, NGC 1596, is very similar 
to the antitruncations presented by lErwin et"al] |2005). The outer 
profile is distinctly shallower in slope than the main inner disk pro- 
file. Although it should be classified as Type III in this sense, the 
profile is not well described by our simple broken exponential fit- 
ting function. The whole disk profile (beyond R^, ± 25") appears 
to be almost continuously upbending and one may argue for an 
extended, exponential break region between breaks a t R « 70" 
and Rx 130". Nevertheless. lErwin et all d2005il2006T) show some 
early-type, face-on counterparts with very similar profile shapes. 

4.1.2 Face-on vs. Edge-on: quantitative comparison. 

As discussed in Section|2]our sample is not unbiased. However, the 
unusually high frequency of galaxies with Type II profiles suggests 
that we indeed have selected against the two other types. For the 
nine Type II galaxies in our sample, we find that the break radius is 
in the range between 1.2— 3.3/ii n • To compare this with the results 
of the face-on sample, we have to account for the effect of dust on 
the inner scalelength. In Section |3T6l we have shown that in the case 
of dusty edge-on galaxies the resulting scalelengths are typically 
20% overestimated (see Fig. |3). Factoring this correction^ in one 
obtains a break-ra dius in the range of 1.4 — 4.0fcj n , which matches 
the range given bv lPohlen & Truiillol J2003) . 

The surface brightness at the break radius ranges from 21.9 
to 26.3 R-mag/rj" , which also match e s the full range for face- 
on galaxies given bv lPohlen & Truiillol < f20o3) . Here, the surface- 
brightness for each galaxy is corrected for galactic extinction and, 
in the case of the V-band images, transformed to the R-band assum- 
ing {V-R) =0.7. 

4.2 Vertical distribution 

4.2.1 The inner vertical structure 

First we concentrate on the inner disk region, i.e., the region in be- 
tween the bulge/bar component and the location of the break. For 
almost all of the late-type galaxies (Sbc and later, for details see 
Appendix [At we find an increase in the measured scalelength with 
the distance from the plane. The increase, at 3 — 5zo, is typically of 

2 Without the small correction for the different wavelengths used (R and V 
band here compared to r'-band for the face-on sample). 
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33.1 
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2.0 
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R 


II 


15 


100 
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27.9 
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24.0 


0.07 
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R 


11 


10 


59 
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5.7 
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2.7 


25.3 
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V 


II 


20 
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40.4 


9.7 


2.9 


4.2 


24.4 


0.27 


FGC2339 


R 


II 


37 


66 


57.3 


17.6 


5.9 


3.3 


3.0 


26.4 


0.10 
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V 


II 





92 


75.6 


34.0 


9.6 


2.2 


3.6 


24.6 


0.13 



Table 3. Exponential disk parameters for the face-on-equivalent profiles. 

(1) Principal name, (2) filter, (3) profile type (I: no break, II: downbending break, III: upbending break), (4,5) inner and outer fitting boundaries, (6.7.S) break radius, 
inner, and outer scalelength in units of arcsec, (9) break radius in relation to the inner scalelength, (10) inner scalelength in relation to the outer scalelength, (11) 
the surface bright ness at the brea k radius (estimated at the crossing point of the two exponential fits), (12) galactic extinction for the appropriate band from 
NED according to Schlegel et alj <1998l) . 



the order of about 20%, relative to the major axis scalelength (see 
upper panel of Fig. [6). Since we know that the omnipresence of 
dust extinction will mostly affect the major axis profiles, as evident 
from the systematic flattening one sees for our artificial galaxies 
(Fig.O, our measurements of the increase in scalelength in the op- 
tical, starting at the major axis, are only lower limits. Consequently, 
the major axis slope for four galaxies is indeed observed to be flat- 
ter than that obtained from adjacent z cuts (cf. Fig.[6]l. 

It is interesting to note that this type of flattening towards the 
midplane, which is for R-band images well explained by the dust 
effect, is also present at the level of approximately 15% in our 
K'-band image (NGC 5290). If this is not an intrinsic feature of 
NGC 5290, one could argue that the dust may still influence the 
scalelength determination on the major axis even in the NIR band. 
On the other hand, this seems unlikely since the amount of dust in- 
ferred from the K'-band flattening should have a more severe effect 
on the R-band profile (see Appendix [At. 

In addition, tests with NGC 5981 (see Appendix [At show that 
using the mean quadrant instead of the dustfree side, also results 
in a systematic underestimation for the measured increase in scale- 
length with distance from the plane. 



4.2.2 The outer vertical structure 

Here we focus on the vertical structure of the outer disk, i.e., the 
region beyond the break radius. For all of the late-type galaxies, 
namely, Sbc and later, a significant increase of the measured scale- 
length with distance from the plane is observed (see lower panel of 
Fig. |6). Whereas for the inner disk the increase in the slope is by 
a factor of 1.1-1.4, for the outer disk we find a factor of 1.4-2.4. 
This indicates an extreme flattening of the profile towards larger 
distances from the plane. Again, the dust, well known for flatten- 
ing surface brightness profiles, cannot be responsible here since it 
is concentrated in the plane. This could be regarded as a weaken- 
ing of the truncation with increasing distance from the plane (see 
Fig.0. 



4.2.3 The break radius as a function of z 

The position of the break radius as a function of vertical distance 
from the plane (z) is another constraint that any theory that might 
explain the truncation phenomenon should satisfy. An inspection 
of the location of the break radius from the deprojected profiles 
(Fig.© shows that this location does not change significantly as a 
function of z in the first two to three z-slices. For the more distant 
z-slices the slopes of the inner and outer regions become closer, 
due the more substantial flattening in the outer disk, which makes 
it harder to confidently determine the location of the break radius. 
A clear demonstration of this is when the inner and outer slopes 
are very close, in such a case the presence of any genuine small 
scale structure in the image renders the location of the crossing 
point very uncertain. This is even more severe when the broken 
exponential model does not provide a good description of the data, 
e.g., when the break region is clearly extended (ESO 380-019 or 
ESO404-018). 

Resorting to the NIR image or concentrating on the smooth 
SO galaxies in our sample, in which the substructure (star-forming 
region or dust patches) are less prominent, does not improve the 
matter because for these galaxies, as with the other, the flattening 
of the outer disk dominates the uncertainty in the location of the 
crossing point. 

To summarise, based on the first 2 to 3 vertical slices, the cur- 
rent data suggest that the location of the break radius for all our 
galaxies is consistent with being constant. 



5 DISCUSSION 
Profile classification 

In our pilot sample of edge-on galaxies we do indeed find the 
sa me three main types of su rface brightness profiles as cla ssified 
by lErwin et all d2005l 120061) and IPohlen & Truiillol d2006h from 
their analysis of face-on galaxy images. The broken exponential pa- 
rameters for the face-on-equivalent view obtained from our edge- 
on galaxies are well consistent with those found by the analysis 
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Figure 5. Deprojected radial profiles of the sample galaxies. 

The surface brightness profiles of each galaxy are shown in separate panels. The different line types are explained in Fig. [4] 



of, the more complete, genuine face-on galaxy samples. Only for 
FGC 2339 the break seems to be at a slightly lower than usual sur- 
face brightness, albeit at a typical distance from the center. How- 
ever, this is also the galaxy with th e value of the lowest ce ntral 
surface brightness. Hence, similar to Pohlen & Trujillg] j2006t) . we 
still cannot disentangle which of the two statements is correct: the 
surface brightness at the break radius is constant, or it is coupled to 
the central surface brightness. 

We do not find a difference in the shape of the truncation com- 
paring the two profiles for NGC 5920 obtained from the optical and 
the NIR image. Both show clearly a broken exponential structure in 



the regi on between 30 ^ R SS> 1 10" . This is consistent with the find- 
ings of iHunter et ail (2006) for the profile shapes of several irreg- 
ular galaxies using 3.6/i and 4.5/i images from the Spitzer Space 
Telescope. Irregular galaxies with double-exponential optical light 
profiles have the same in the NIR. So baring an exception, NIR 
observation of truncations do not seem t o differ significantly from 
observation in the optical as proposed bv lFlorido et al.l d2006T) . 

Pohlen &Truiilld ( l2006t) have shown tha t Type I galaxies ar e 
rare in late-type galaxies and, according to lErwin et all d2006h . 
more frequent in SO galaxies. Accordingly, our only Type I can- 
didate (NGC 3390) exhibits an SO-like outer disk structure, despite 
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being classified as Sb. It is interesting to note that in this case the 
profiles of the individual cuts, close to the major axis but uninflu- 
enced by the dustlane, show some indication of being better de- 
scribed by a broken exponential structure. Therefore, the vertical 
integration in this case may obscure the existence of a break close 
to the plane. This is clearly caused by the lack of an apparent break 
in the vertical slices at larger vertical height above the midplane. 
This might explain the high frequency of not-truncated disk s seen 
in the early-type face-on galaxy sample of Erw in et aT] J2006h . 

Despite the discussion above about the general agreement be- 
tween the edge-on and face-on classifications, an edge-on coun- 
terpart of the not-truncated, Type I , late-type galaxies (e.g., the 



face-on NGC2776 bv lPohlen & Truiilloll2006l) . is yet to be found. 
Given the expected frequency of such galaxies (~ 10%) and the 
size of the parent sample from which we select, we are not able to 
explain the absence of such galaxies. 

If late-type Type III galaxies might be a ssociated with in- 
teracti ng or merging systems, as proposed by IPohlen & T ruiillo 
( 2006), the detection of only one Type III galaxy is in full agree- 
ment, as we excluded those systems from the original, late-type 
dominated, sample. In addition, the fact that our Type III galaxy 
is an early-type one is consistent with the observed increase in 
frequ ency of Type III profiles for early-type galaxies dErwin et al.l 
2006). To our knowledge, all edge-on Type III galaxies reported in 
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Figure 6. Measured scalelength with increasing distance from the plane. 
Upper panel: The open circles are the measured inner scalelengths (h; n ) 
normalised to the inner scalelengths measured on the major axis /iin ma 
for each profile, versus the distance (z) from the plane normalised to the 
scaleheight (h z ) given in Appendix [A] The filled squares mark the values 
for the outermost profiles. All points for an individual galaxy are connected 
(solid lines for Type II, dotted-dashed line for Type I, and dashed line for 
Type III) and the name of the galaxy is added. 

Lower panel: The open circles are the measured outer scalelengths (h ou t ) 
normalised to the outer scalelengths measured on the major axis ftout ma 
for each profile versus the distance (z) from the plane normalised to the 
scaleheight (h z ). 



the literature so far are ear ly-types. These galaxies are NGC 4762 
studied by several gr oups jTsikoud3ll98d : [Wakamatsu & Hamabe l 
ll984l:IWozniakll994h. NGC 4452 shown bvlrlamabe & WakamatstJ 



( 1989) and NGC 4696C0 fromlPohlerJ j200lF ) 
For face-on galaxies. lErwin et ali i 2005) 



differentiate between 
two different types of Type III profiles. One where the excess light 
is probably part of the disk which show a relatively sharp transi- 
tion and both parts are well described with exponentials (Type III- 
d). And another one where the excess light is associated with the 
halo/spheroid or bulge (Type III-s). The latter is characterised by 
a small inner region fitted with an exponential and a continuously 
upbending outer profile and a changing ellipticity with radius. In 
the case of NGC 1596 we would argue for a Type III-s classifica- 
tion where the outer structure is associated more with an envelope 

3 It is worth noticing, tha t NGC4696C is classified as an Sb galaxy 
jde Vaucouleurs et alJI 199 lh . However, deep imaging, obtained bv lPohlen] 
12001), clearly shows an outer SO-like envelope. 
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Figure 7. Weakening of truncation with increasing distance from the plane. 
The open circles are the measured ratios of inner to outer scalelength 
(hi n /h ou t ) for each Type II profile with increasing distance (z) from the 
plane normalised to the scaleheight (h z ). The filled squares mark the values 
for the outermost profiles. All points for an individual galaxy are connected 
(solid lines) and the name of the galaxy is added. 



like structure. This is supported by the difficulty one faces in fit- 
ting NGC 1596 wit h a consistent thin/thick disk combination (cf. 
iPohlen et alj|2004bf) . Any bulge component would have to have an 
unlikely, extremely high, axis ratio (cf. Appendix [Aj. In this case, 
the observed outer envelope is ther efore much better e xplained by 
interaction as recently supported by Chun g et al. (2006). 

Inner Disk 

Almost all of the analysed profiles exhibit an increase of the ra- 
dial scalelength with increasing distance from the plane. Notwith- 
standing the lack of completeness in our sample, this finding 
co uld support the notion of ub iquitous thick disks as proposed 
bv lDalcanton & Bernstein! {2002), where these thick disks have to 
have a systematically larger scalelength, compared to the thin disk. 
However, in the absence of a consistent two disk fit to the data at all 
radii the thick disk model is not the only option. Instead one can in- 
terpret the increase in the radial scalelength merely as a result either 
of a radially flattened component above the inner disk or a flatten- 
ing of the inner disk itself. For example, in the case of the late-type 
(Scd) galaxy FGC 2339, the whole inner disk region is dominated 
by a ring- or bar-like structure in the midplane which also alters the 
shape of the disk above/below the plane over the entire measured 
vertical extent. 

Outer Disk 

In order to examine what could lead to the observed flattening of the 
outer disk profile, we invoke several 3 dimensional (3D) toy mod- 
els. These 3D models are generated by combining separate compo- 
nents, e.g., a truncated thin disk plus a truncated thick disk. The re- 
sultant structures are qualitatively compared to those reconstructed 
from the data. It is important to emphasise here that no attempt 
is made to fit the data to these models, rather the intention is to 
explore whether these toy models could generically reproduce the 
increased flattening seen in the data. The 3D profiles of the models 
are shown in Fig. [8] 

To highlight the flattening in these models we fit broken expo- 
nential profiles with the same fitting routine used for the real galax- 
ies (cf. Section |33t . For example the upper-left panel (1) of Fig. [8] 
shows a single broken exponential disk model. This type of models 
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Figure 8. Toy models to explain observed profile structure. 

Deprojected radial surface brightness profiles of six model galaxies: single disk (panel 1), truncated thin plus not-truncated thick disk (panel 2), and two 
models with flaring (panels 3-4). The figure has the same setup as Fig. [5] 



does not show any flattening in the reconstructed scalelengths and, 
therefore, is not favoured by our data. 

Could we be seeing de V aucouleurs R 1 ^ bulges? 
iBahcall & Kvlafij dl985h have shown that a thin disk plus a R 1/4 
bulge can mimic the shape of the isophotes of a thin disk plus thick 
disk model (without bulge) down to a surface brightness level of 
about hr ~ 25.5mag/ □" . Furthermore, this combination can even 
produce some flattening of the radial profiles with increasing verti- 
cal height beyond the break radius as seen in our models. 

Having said that however, there are strong arguments against 
explaining the observed flattening with a spherical power law bulge 
component. F i rst, in a sample of 19 early-type galaxies (SO-Sbc) 
iBalcells et aT] J2003b find that none of the bulges have a i? 1/4 be- 
haviour. Most of their bulges are better fit with exponentials which 
are not expected to dominate the radial profiles at large distances. 
Therefore, it is statistically unlikely to attribute the flattening in the 
outer disk in most of our late-type edge-on galaxies to a spherical 
bulge component. Second, the observed shape of the outer contours 
in all the cases studied here appear too much disk-like to be caused 
by a spheroid (see Fig. IB II and Fig. [9). Third, it is almost impos- 
sible to create a _R 1//4 bulge which will leave the high S/N major 
axis profile purely exponential down to the last measured point, but 
will dominate the radial profile at larger vertical heights above the 
plane, and doing this all without being clearly visible in a contour 
map. 

In conclusion, since there is no clear sign of a de Vaucouleurs 
i? 1 / 4 in any of the observed images whereas the outer disk flatten- 
ing is found in the majority of galaxies in the sample, the answer to 
the question posed here is, most probably: no. 



Does a thin plus thick disk combination produce the increased flat- 
tening beyond the break? 

If one assumes that the thick disk scalelength is in general larger 
than the thin disk scalelength, the increase in scalelength with z- 
height in the inner disk will be always present. However, there are 
various possibilities for producing the extended flattening in the 
outer disk beyond the break with thin plus thick disk models. The 
exact shape of the flattening critically depends on the, unknown, 
outer disk properties of the thick disk. 

If, for example, we assume the thick disk to be not-truncated, 
we are still able to produce some flattening beyond the break (see 
panel (2) of Fig. [8J. However, depending on the surface bright- 
ness limit the not-truncated thick disk will become visible at large 
radii. This could be an explanation for some of the very few 
Type II + III mixed clas sifications found for late-type galaxies by 
IPohlen & Truiillol J2006h . These galaxies show two breaks in their 
face-on surface brightness profile. The first is downbending in the 
inner region, while the other is upbending in the very outer re- 
gion. However, none of our edge-on sample galaxies shows this 
behaviour. 

By adding a truncation to the thick disk one can produce a 
wide variety of broken exponential structures depending, for ex- 
ample, on the ratio of inner to outer disk scalelength (hi n /h out )• 
Models with unequal thin and thick disk truncation radii provide 
even more choice to fit the data. 

In conclusion, the answer to the question posed here is: yes. 
However, the degree of flattening will depend on the details of the 
3D model. 
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Figure 9. Upper panel: Contour map of a radially broken exponential, ver- 
tically exponential thin disk superimposed with a standard de Vaucouleurs 
R 1 / 4 bulge. Lower panel: Three combination of a the same thin disk and a 
thick disk, where the thick disk scaleheight grows (flares) between the break 
radius (7?br ) an d three times with an increasing factor (solid, dashed, 
dotted lines). 



Could the flattening be explained by intrinsic flaring_\ of the stellar 
disk? 

To produce an increased flattening in the outer-disk one can invoke 
models with i ntrinsic flaring - similar to those used for the Milky 
Way (see e. dGvuket aiT l999) - that starts beyond the break ra- 
dius. In the edge-on view, the contour map of the observed galaxy 
provides a strong constraint on the degree of flaring. The lower- 
panel of Fig. [9] shows three different radially broken exponential 
single disk models with flaring factorjf] of F = 2.5, 5, and 10. If 
the flaring factor is chosen to be too large the resulting contour map 
will blow up in the outer parts, a behaviour that is neither observed 
in any of our sample galaxies nor in the much larger sample of 
In general, it seems that a strong flaring of the stellar 
disk is rarely observed. In addition, in the few reported cases it ap- 
pears to be rather associated with Type III profiles, e.g., NGC 476 2 
i lWakamatsu & Hamabdl 19841 ; iFreeman & Bland-Hawthorr3l2002h . 
and not with Type II as discussed here. 

As shown in the panel (3) of Fig. [8] flaring of a single disk 
beyond the break radius does not always produce the right amount 
of flattening if compared to our observations. A more promising 



4 With intrinsic flaring we do not mean t he sometimes observed fact of an 
increase in scaleheight with radius (e.g. lde G riis & Peletiei 1997), which is 
most probably related to a combination of thin and thick disk with different 
scalelengths, but the intrinsic change of the scaleheight for a single stellar 
disk component. 

5 For the Milky Way iGvuk et~aT] jl999l) have models with F = 5 or 10 
between the solar radius and 25 kpc. 



attempt to apply flaring is obtained by allowing only the thick disk 
component to flare beyond the break radius (see panel (4) of Fig. [8}. 
This approach could help in cases where the inner disk indicates 
only a very small value for the ratio of thin to thick disk scalelength. 
The flaring is then needed to reproduce the observed significant dif- 
ference of the outer scalelength on the major axis compared to high 
above the plane. A flaring of the thick disk is also p roposed by sev- 
eral authors (e.g. Isiegel et al]|2002k |Pu et al.ll2006l) for our Milky 
Way from star-count studies. This flaring would account for the ob- 
served variable scaleheight with different direction of the analysed 
fields. 

It is clear that the data can in principle discriminate between 
possible 3D models, as each will produce different features in the 
reconstructed profiles. However, we make no such attempt in this 
paper except to point out the different possibilities. 

Break Radius 

As we have shown earlier, the location of the break-radius in the 
current sample is found to be consistent with being height indepen- 
dent. In light of the enhanced flattening of the outer disk relative to 
the inner disk, the determination of the break radius at large vertical 
cuts is uncertain and more accurate modelling could reveal signif- 
icant variation in its location towards the larger distance from the 
major axis. 

Nevertheless, we would like to point out that if such a finding 
is supported by more data, it could have far reaching implications 
on theories that attempt to understand how the thin and thick disk 
are assembled. As such a theory should explain how potentially 
distinct components of the disk "coordinate" the location of their 
break radii. 



6 CONCLUSIONS 

This paper has presented a new algorithm for reconstructing the in- 
trinsic 3D structure of galactic disks seen edge-on. The algorithm 
assumes axial-symmetry and utilises the so-called Fourier Slice 
Theorem which in the case of edge-on galaxies is equivalent to an 
inverse Abel-transform. The algorithm is applied to a pilot set of 12 
edge-on galaxy images, drawn from a parent sample of 72 galaxies, 
to obtain face-on-equivalent radial surface brightness profiles, and 
reconstruct the radial profiles along several vertical cuts, including 
the major axis. 

Comparing the results obtained from the face-on-equivalent 
profiles with those measured from "real" face-on samples we reach 
several conclusion. The first, and probably most important, con- 
clusion is that we find in the current edge-on sample at least one 
example for each of the three basic profile types found in face- 
on samples: no break, upbending break, an d downbending break 
dErwin et ID 12001 120061 ; iPohlen & TraUflcl l2006h . Despi te this 
agree ment, our morphologically selected parent sample |Pohlen 
l200lh . choosing undisturbed and prototypical disk galaxies, clearly 
favours late-type galaxies with class Type II profiles, i.e., down- 
bending breaks. This might be due to a connection between certain 
profile types and distinct vertical structures. For example, if late- 
type Type III galaxies are connected to interaction and/or merging 
processes, they would probably be excluded from edge-on sam- 
ples. In other words, they would exhibit distortions, lopsidedness, 
or warping, which are all prominent in the edge-on view. 

Although our parent sample is not unbiased in this sense, our 
results do show that the parameters for the downbending breaks of 
the face-on-equivalent profiles obtained from our edge-on galax- 
ies are consistent with those seen in face-on samples. This implies 
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that truncated profiles of galaxies are independent of geometry and 
there is no distinction in this regard between edge-on and face-on 
studies. 

In this paper we have been able to study, for the first time, 
the deprojected radial profiles as a function of vertical height above 
the plane, which is only possible to explore with edge-on galaxies. 
Hence, this approach provides a unique opportunity to give addi- 
tional constraints on theories explaining the origin of the various 
types of observed profiles in the context of disk formation and evo- 
lution. 

In this context, we find a generic flattening of the observed 
profiles in the inner disk region with increasing distance from the 
plane, which cannot be explained by dust extinction. This implies 
a more complex vertical structure in the inner region than a single 
disk. It is possible to explain this flattening by introducing a thick 
disk with systematically larger scalelength compared to that of its 
host thin disk. 

In the outer disk, beyond the break radius, we find an enhanced 
flattening of the deprojected profiles, which again sets strong con- 
straints on theories that attempt to explain the observed broken ex- 
ponential structures. Using some toy models we conclude that thin 
plus thick disk combinations are, in principle, able to explain the 
observed structures. However, the details of this combination de- 
pend on the exact shape of the outer region of the thick disk. For 
example, a thick disk without truncation will exhibit less flattening 
in the outer disk than a thick disk with truncation. So far our data 
does not allow us to decide if thick disks are in general truncated 
or not. It remains unclear whether one can fine tune a simple model 
with broken exponential thin plus thick disks in order to explain the 
observed brightness distribution. 

In addition, we find that the location of the break radius does 
not change significantly as a function of z in the first two to three 
z-slices. Unfortunately the flattening of the radial profile beyond 
the first few vertical cuts makes it hard to determine the location of 
the break radius. Although we do not know yet whether thick disks 
are indeed truncated, we have shown that the vertical structure har- 
bours clues to understand the origin of the different radial profiles 
and should be pursued. For example, luminous, not-truncated thick 
disk components could be responsible for the high frequency of 
not-truncated lenticular galaxies in face-on surveys, while having a 
broken exponential structure close to the midplane of the galaxy. 

In conclusion, observations of the vertical structure set strong 
constraints on modelling/simulating disk galaxies and the interpre- 
tation of their various components. In the future, a larger, well se- 
lected sample of edge-on galaxies will be analysed in order to shed 
light on some of the issues raised here. 
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APPENDIX A: NOTES ON INDIVIDUAL GALAXIES 
ESO 380-019: 

Owing to a not well defined thick patchy dust lane the inclination 
is uncertain (clearly not 90°, but closer to 85°). The extent and 
size of the bulge component is also uncertain. The major axis ra- 
dial profile shows a broken exponential structure with a break at 
R w ±(95-100)" (fiv « 23.8mag/ Q " ). The galaxy is moder- 
ately asymmetric. In the cont our map the SE (left)-side is obviously 
shorter. The single disk fit bv lPohlenl ( l200ll) yields a vertical scale- 
height of h z — 3.5", a radial scalelength of h = 22.8", and an 
inclination of 84.0°. 

The face-on-equivalent profile of ESO 380-019 is well described 
as a Type II with an inner shallow exponential (of hi n = 27.9" = 
2.4/i ou t ) and an outer steeper exponential (h out = 11.5") with a 
break at a radial distance of i?br ~ 88" = 3.1/iin and a surface 
brightness of /ib r = 22.9V-mag/[]" . 

We find an increase in scalelength of m 25% between the major 
axis and about 5 times the vertical scaleheight (at 18") in the in- 
ner disk region. Since the galaxy is rather far from edge-on, the 
extended dustlane (with maximum extinction at « 7") could still 
play a role (cf. Fig. [3]>. The scalelength beyond the break in the 
outer disk increases by a factor of 2.5 between and 5.1-/i z . 

ESO 404-018: 

This galaxy shows hardly any bulge component . The disk struct ure 
is quite knotty and appears bended like a wave dWeinberdl200lh in 
addition to being slightly u-shaped at the edges. The radial surface 
brightness profile shows quite clearly the deviation from a simple 
not-truncated exponential model on both sides with a sharp break 
on the NE(left)-side, and a smoother one on the SW(right)-side 
starting at 7? ~ ±70 " and is well fitted by a broken exponential. 
The single disk fit by Pohlerj J200lh yields a vertical scaleheight of 
hz — 2.8", a radial scalelength of h = 31.3", and an inclination of 
90.0°. 

The face-on-equivalent profile of ESO 404-0 18 is well described 
as a Type II with an inner shallow exponential (of h ln = 34.0" = 
3.6/iout ) and an outer steeper exponential (h out = 9.6") with a 
break at a radial distance of i?br ~ 76" = 2.2h ln and a surface 
brightness of /ibr = 24.6V-mag/n" . 

We find an increase in scalelength of « 12% between the major 
axis and about 4 times the mean vertical scaleheight (at 11") in the 
inner disk region. Although this galaxy is almost perfectly edge-on 
we still find the steepest profile on the major axis, so the dust is not 
playing a major role here. The shape of the profile in the outer disk 
is not well described with a single exponential. The break regions 
(starting at w 75") seems to be extended and almost exponential 
itself. This deviation from our simple broken exponential model 
strongly affects the measurement of the scalelengths beyond the 
break. 

FGC 2339: 

This galaxy shows on both sides slightly off-centered ring-, or bar- 
like knots at (R,z) = (±24", ±1") and (R,z) = (-26", ±1"). 
The underlying disk is quite asymmetric on both sides. The ma- 
jor axis profile on the N(left)-side is satisfactorily fitted with a 
broken exponential showing a break at R ~ ±(55 — 60)" and 
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~ 24.0mag/[]" . On the S(right)-side the profile looks more 
like a sharply-truncated profile. The single disk fit bv |Pohlenll200lh 
yields a vertical scaleheight of h z — 1.8", a radial scalelength of 
h= 19.9", and an inclination of 88.5° . 

The face-on-equivalent profile of FGC 2339, beyond 35" can be 
described as a Type II with an inner shallow exponential (of hi n = 
17.6" = 3.0/iout ) and an outer steeper exponential (h out = 5.9") 
with a break at a radial distance of Rh T ~ 57" = 3.3h[ n and a 
surface brightness of /ib r = 26.4R-mag/ □" . 
For this galaxy the whole inner region R Ss 35" is dominated 
by a ring or barlike structure which introduces an extreme flat- 
tening of the profiles a t 2 vertical scaleheights and above (cf. 
Ivan den Berg et ai]|2006h . The small inner disk region in the range 
of 35" R 5s 55" does not show a clear systematic flattening to- 
wards higher z-cuts. The major axis profile is overall the flattest, 
so probably heavily affected by the midplane dust. Nevertheless, 
there is a 19% increase in scalelength between 2.3 and 3.4-/i z . The 
scalelength beyond the break in the outer disk increases by a factor 
of 2.3 between — 2.3- h z . 

IC4871: (alternative name: IC4872) 

This galaxy is overall very patchy with many distinct H II regions. 
On the S(right)-side, along the major axis at 7? = +100" there 
is an unusually warped patch of star formation visible. The major 
axis radial profile can be well described with a broken exponen- 
tial with a very steep outer decline (dropping by « 2 mag within 
« 10") on both sides (with slightly different outer slopes), start- 
ing at R ~ ± 115"(/iy ~ 23.5mag/n"). The single disk fit by 
Pohlen J200lh yields a vertical scaleheight of h z — 3.8", a radial 
scalelength of h = 42.3" , and an inclination of 87.5°. 
The face-on-equivalent profile of IC4871 is well described as a 
Type II with an inner shallow exponential (of hi n — 40.4" = 
4.2ftout ) and an very steep outer exponential (h out =9.7") with a 
break at a radial distance of R^ T ~ 117" = 2.9h m and a surface 
brightness of /Xb r = 24.4V-mag/n" . 

We find an increase in scalelength of w 13% between the major 
axis and about 2.5 times the mean vertical scaleheight (at 8") in 
the inner region. The scalelength beyond the break in the outer disk 
increases by a factor of 1.7 between — 2.5- h z . 

NGC522: 

Our image has some residual structure in the Hatfield which reaches 
the disk on the SW(right)-side at (R,z) = (+10", -15") and at 
(R,z) — (+60",— 30") causing the unusually shaped contours. 
There is a galaxy cluster in the background with an edge-on spi- 
ral galaxy superimposed at (R,z) — (—13",— 17") and another 
elliptical at (R, z) = (—60", —20"). The major axis radial pro- 
file shows a clear broken exponential behaviour with a very steep 
outer decline (dropping by~ 5 mag within m 25") on both sides, 
starting at R « ±60" (fin ~ 21.5mag/[]" ). Even the edge-on 
profiles show clearly a decrease in the outer slope with increasing 
z (see below). The galaxy is slightly asymmetric with the N E(left)- 
side sh orter than the SW(right)-side. The single disk fit bv lPohlenl 
yields a vertical scaleheight of h z = 4.3", a radial scale- 
length of 7i = 32.0", and an inclination of 87.5°. 
The face-on-equivalent profile of NGC522 is well described as 
a Type II with an inner shallow exponential (of h- m = 52.7" = 
6.1ftout ) and an very steep outer exponential (h out = 8.7") with 
an 'early' break at a radial distance of R^ r ~ 63" = 1.2hi n and a 
surface brightness of /ibr =22.8R-mag/n" . 
The major axis profile is overall the flattest, so probably heavily 
affected by the midplane dust. We find only a marginal increase in 
scalelength of m 5% between 0.9 and 3.1 times the mean vertical 



scaleheight (between 4 — 13") in the inner region. The scalelength 
beyond the break in the outer disk increases by a factor of 1.8 be- 
tween — 3.1- h z . 

NG C 1596: 

See iPohlen et al. | d2004bl) for additional notes. The face-on- 
equivalent profile of NGC 1596 is classified as Type III, i.e., a break 
with an upbending profile beyond. However, the profile is not well 
described by a simple broken exponential. The whole disk profile 
(beyond R^ ± 25") appears to be almost continuously upbending. 
One may argue for an extended, exponential break region between 
breaks at R « 70" and R as 130". Therefore, the exact position 
of the break is unclear. Judging from the individually deprojected 
cuts this structure is driven by the clear upbending breaks on the 
major axis and at z—8" at a radial distance of 110". All cuts at 
larger distance from the plane (z<;17") are consistent with having 
no break. 

We find an increase in scalelength by a factor of 3.6 from the major 
axis (steepest) to a vertical scaleheight of 33" (flatter) in the inner 
region. 

NGC 3390: 

There is a nearby, low surface brightness companion, of unknown 
distance, visible at (R, z) = (+140", —125") and another unusu- 
ally small and round object at R = +85" on the major axis. This 
galaxy exhibits a quite unusual appearance. The inner high sur- 
face brightness part looks like a normal edge-on Sb galaxy with 
a well defined dust-lane. However, the outer low surface brightness 
part shows that the galaxy is embedded in a huge spheroidal enve- 
lope mu ch more l i ke an SO. The galax y was already classifi e d as a n 
SO-a bv lLaubertsI dl982l) and SO? bvlSandage & Tammannl dl980t). 
whereas de Vaucoule urs et alj dl99 1 ) and Prugniel & Heraudeau 



dl998l) quote Sb or Sbc. The radial profiles close to the major axis, 
similar to other SO galaxies, decline almost straight into the noise 
with only a mild break at R m ±90" in cuts higher above the plane. 
This structure makes it impossible to fit a sharply truncated or bro- 
ken exponential model. Even the edge-on profiles show clearly a 
decrease in the outer slope with increasing z (see below). 
The face-on-equivalent profile of IC 3390 is reasonably well de- 
scribed as a Type I (see discussion in main text, Section |4~TT >. For 
NGC 3390 we find an increase in scalelength by a factor of 3.1 
between the major axis (steepest) and a vertical distance of 26" 
(flatter) in the inner disk region as expected from the SO envelope . 
A scaleheight of 5.3", measured in the I-band bv ldeGriisI d 19981) . 
is used in Fig.[6]and Fig.[7] 

N GC 4179: 

See lPohlen et"afl (2004b) for additional notes. We find that the face- 
on-equivalent profile of NGC 4179 is well described as a Type II 
with an inner shallow exponential (of h- m = 30.9" = 1.4/i ou t ) and 
an outer slightly steeper exponential (h out = 22.6") with a break 
at a radial distance of R^ T w 98" = 3.2ftj n and at a surface bright- 
ness of /ibr =23.3V-mag/o" . 

In the inner disk we find an increase in scalelength of 36% between 
the major axis and 4.5 times t he thin disk scaleheight (6.1") re- 
ported by IPohlen et all (2004b ) . Overall the scalelength increases 
by a factor of 1.8 to the highest measured c ut at z = 45" . This is 
very similar to the factor of 1.9 reported by IPohlen et alj (2004b) 
for the ratio of thick to thick disk scalelength, obtained with a com- 
pletely different method (3D model fitting). Their thin disk scale- 
length of 26.2" is almost identical to our scalelength of 26.7" at 
the major axis. Our scalelength of 46.8" for the vertical cut with 
the largest distance from the plane (at z — 45") is » 9% smaller 



20 M. Pohlen, S. Zaroubi, R.F. Peletier, & R. -J. Dettmar 



compared to their thick disk scalelength of 51.0", which we may 
reach for cuts even higher above. 

NGC5290: 

The inner disk is embedded in an outer, clearly thicker disk com- 
ponent or in a low surface brightness extension of the bulge. The 
galaxy exhibits a thick and patchy dustlane, and the two sides are 
rather asymmetric. In the contour map the E(left)-side seems to be 
thicker and shorter, whereas the W(right)-side is longer and tapers 
off. This asymmetry is reflected in the major axis radial profile. 
When fitted with a broken exponential structure, the slope of the 
outer region on the E(left)-side is much steeper compared to the 
W(right)-side. The latter seems to be almost ~ 15" 'longer'. The 
single disk fit to the R-band image bv lPohlenl j200lh yields a verti- 
cal scaleheight of h z = 7.5", a radial scalelength of h = 25.4", and 
an inclination of 88.5°. 

The face-on-equivalent profile of NGC 5290 is in both bands (R 
and K') well described as a Type II with an inner shallow exponen- 
tial (of h in = 33.1" = 2.0/w in R and h in = 31.2" = 2.3h out 
in K') and an outer steeper exponential (h out = 16.3" in R and 
h ou t = 13.6" in K'). The break radius is at a similar radial distance 
(R hr S3 80" ee 2Ah- m in R, and R b[ S3 83" = 2.7 h in in K') and at 
a surface brightness of /Xbr = 23.7R-mag/n" and /ibr = 21.0K'- 
mag/n" respectively. 

In the K'-band, the scalelength of the major axis in the inner disk 
region is the flattest compared to all other vertical cuts above (see 
Fig. O. For the R-band it is almost as fiat as the scalelength at a 
distance of 4 times the vertical scaleheight. We observe a system- 
atic decrease in scalelength in both bands up to a vertical distance 
of z S3 15", after which the scalelength starts to rise again for the 
deeper R-band image out to the highest vertical cut at 30" . If this is 
not an intrinsic feature of NGC 5290 one could argue that the dust 
may still influence the scalelength determination even in the NIR 
band. However, quantitatively this seems very unlikely since the 
amount of dust inferred from the K'-band flattening should have 
a more severe effect on the R-band profile. In the deeper optical 
image the scalelength increases by 26% between 2 and 4 vertical 
scaleheight. In the R-band the scalelength beyond the break in the 
outer disk increases by a factor of 1.5 between — 2.6-/i z . 

NGC 5981: 

There is a possible dwarf elliptical companion visible at (R, z) = 
(—10", +30"), and another companion superimposed on the disk 
at (R, z) = (+77", —5"). For both we have no velocity informa- 
tion, but there is another large elliptical nearby (NGC 5982, 6.2' 
away with At; S3 1300 km s" 1 ). NGC 5981 is slightly disturbed and 
asymmetric. The inner disk is S3 1.5° tilted against the outer disk 
and the bulge component is furthermore slightly tilted against the 
inner disk. This asymmetry is reflected in the radial profiles. The 
outer slope of a broken exponential fit is significantly shallower on 
the NW(right)-side, where also the contour map is more extended. 
The single disk fit by Pohler] j200lh yields a vertical scaleheight of 
h z — 4.2", a radial scalelength of h = 24.1", and an inclination of 
86.5°. 

The face-on-equivalent profile of NGC 5981 is well described as 
a Type II with an inner shallow exponential (of h- m = 27.9" = 
2.0/iout ) and an outer steeper exponential (h out = 13.8") with a 
break at a radial distance of R br S3 67" EE 2Ah[ n and at a surface 
brightness of /ib r = 24.0R-mag/n" . 

We find an increase in scalelength of S3 10% between the major 
axis and about 3.6 times the mean vertical scaleheight (at 15") in 
the inner region. Using only the two quadrants from the dust free 
side for averaging (in opposite to all four) this increase in scale- 



length for the inner disk grows to S3 25%. In this case - as ex- 
pected from the model fits (cf. Fij^ _ the scalelength on the major 
axis is less affected by the dust, so smaller, whereas the scalelength 
at 3.6 vertical scaleheights is almost identical (29.7" for the mean 
quadrant versus 29.2" for the dust free part). The scalelength be- 
yond the break in the outer disk increases by a factor of 1 .4 between 
— 3.6- h z in our mean quadrant. Owing to the described asym- 
metries in the disk on both sides this galaxy is well suited to study 
their effects on the deprojected profiles. Fitting only a mean quad- 
rant obtained from the NW(right)- or SE(left)-side of the galaxy 
separately we find differences of 5% — 10% for the scalelengths of 
the inner region. In the outer disk (beyond the break) the differences 
are clearly larger, e.g., 40% on the major axis, due to the obvious 
asymmetry. 

UGC 10459: 

This galaxy has only a very small bulge component. There is 
a possible background galaxy, with unknown distance, visible at 
(R,z) = (+33", —9"). The major axis profile shows an obvious 
broken exponential structure with a break radius at R S3 ±45" 
(/LtR ss 23.0 mag/ □" ) separa ting a shallow inner and a steeper outer 
slope. The single disk fit bv lPohlenl l200l) yields a vertical scale- 
height of h z — 2.0", a radial scalelength of h — 13.4", and an 
inclination of 87.0°. 

The face-on-equivalent profile of UGC 10459 is well described as 
a Type II with an inner shallow exponential (of h[ n = 15.6" = 
2.7/iout ) and an outer rather steep exponential (h out =5.7") with 
a break at a radial distance of 7?br ~ 44" = 2.8/Hn and at a surface 
brightness of /ibr =25.3R-mag/n" . 

We find an increase in scalelength of S3 40% between the ma- 
jor axis and about 5.6 times the vertical scaleheight (at 11") in 
the inner disk region. The scalelength on the major axis is slightly 
larger compared to scalelength of the cut at 2.1 -h z , which is well 
explained by the influence of the dust. The scalelength beyond 
the break in the outer disk increases by a factor of 2.2 between 
0-3.6-hz ■ 



APPENDIX B: FIGURES 

The following figures show radial surface brightness profiles and 
isophote maps for four galaxies per page (in alphabetical order). 
For the isophote map (lower panels) each galaxy is rotated to the 
major-axis. The magnitude of the outer contour (/iii m , defined by a 
3a criterion of the background) is indicated in each plot. The con- 
secutive contours are equidistant spaced by 0.5 mag. In the plot, 
the contour lines are drawn with increasing smoothing towards the 
outer parts. Note however that for the actual fitting a constant scale 
smoothing is used. For the inner contours out to where the noise 
begins to increase no smoothing was applied. The following 2-3 
contours are smoothed by replacing each pixel by the mean of 3 2 - 
pixels "around" the central pixel. For the outer two contours this 
smoothing is increased to 5 2 -pixels. The upper panels displays the 
major-axis surface brightness profile (top, solid line) and two par- 
allel radial profiles (lower, solid lines) each above and below the 
major axis. The exact vertical positions (z) for the plotted profiles 
are indicated in the upper right corner of the plot. In some c ases the 
best fi t, sharply truncated, single disk model obtained by IPohler] 
(2001) is overplotted (dashed line). 
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Table Al. Exponential disk parameters for all fitted profiles. 

(1) Principal name, (2) filter, (3) profile type (I: no break, II: downbending break, III: upbending break), (4) vertical position z of the deprojected slices ('feo' 
means the face-on-equivalent profile), (5,6) inner and outer fitting boundaries, (7,8,9) break radius, inner, and outer scalelength in units of arcsec, (10) break radius 
in relation to the inner scalelength, (11) inner scalelength in relation to the outer scalelength, (12) the surface brightness at the break radius (estimated at the 
crossing point of the two exponential fits. The full version of this table is available online. 
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Figure Bl. Radial surface brightness profiles and isophote maps for ESO 380-019 V-band, ESO404-018 V-band, FGC2339 R-band, and IC4871 V-band 
(from top left to bottom right). See text for further explanations. 
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Figure Bl. Radial surface brightness profiles and isophote maps for NGC 522 R-band, NGC 1596 R-band, NGC 3390 R-band, and NGC 4179 V-band (from 
top left to bottom right). 
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Figure Bl. Radial surface brightness profiles and isophote maps for NGC5290 R-band, NGC5290 K'-band, NGC5981 R-band, and UGC 10459 R-band 
(from top left to bottom right). 



